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A fully microscopic theory for the spontaneous emission from semiconductors is discussed. The 
theory is evaluated for a quantum-well system and the role of excitonic and unbound electron-hole- 
pair contributions to the emission is analyzed. Simplifying aproximations to the full theory and 
^^ their range of validity are discussed. Numerical solutions are presented for experimentally relevant 

'^'h . situations and it is shown that a detailed analysis of measured spectra requires the knowledge of 

j^ ■ both absorption and emission under identical conditions. 

I. INTRODUCTION 

^ . 

r^ ' In this paper, we present a fully quantum mechanical analysis of spontaneous light emission from a Coulomb 

interacting carrier system in a semiconductor quantum well (QW). We are particularly interested to analyze the 

fl^ ' influence of excitonic and of unbound electron-hole-pair populations on the photo luminescence (PL) emitted at the 

spectral position of the exciton resonances. 

Extensive experimental research in the past few decades has shown that one observes excitonic resonances in the 

"jrt \ PL of a QW for different material systems and for a very broad range of excitation conditions ranging from cryogenic 

(about 4 K) all the way to room temperaturesi^ A PL resonance at the lowest Is-exciton state has also been observed 

'T' ' even when the exciton ionizaton density is exceeded. For example, experiments with GaAs-type QW systems at 4 K 

'^ , lattice temperature show a clear excitonic resonance at the Is position for a density range 10*cm~^ to 5 x lO^^cm"^— ; 

^ ■ here, the lower density is only limited by the sensitivity in detection while the upper value is already beyond the 

O \ ionization density. 

The build-up of excitonic PL is often investigated in a configuration where the semiconductor is excited nonreso- 

nantly using ultrafast laser pulses that generate carriers into the high-energy states of the bands. Thus, all coherences 

vanish almost instantaneously and the direct generation of Is excitons can be excluded. Using this scheme, it has 

^ ■ been observed^ that an excitonic Is PL resonance builds up on a subpicosecond timescale. Furthermore, excitonic 

~V . PL is observed longer than several nanoseconds, which also implies that electron-hole pairs exist in the system for a 

2jL ■ very long time. 

^^ ', In order to explain some of the observe phenomena, Refi^ developed a simple microscopic model for the semicon- 

f^ ' ductor PL showing that also an electron-hole plasma, i.e. a system with nonvanishing electron and hole densities 

\^ \ but no excitons, can produce PL spectra that peak at the exciton resonances. This prediction releases the previously 

^^ I assumed strict relation between the existence of exciton populations and the appearance of excitonic PL such that 

one can explain why excitonic PL is also observed under conditions where the presence of real exciton populations 

seems unreasonable. 

^ ' In order to present a more complete analysis, we show in this paper how one can analyze self-consistently the 

_JL ', electron- hole plasma and exciton-population effects in the QW PL. We show that in order to perform a reliable 

(-H • analysis, one has to evaluate not only the PL but also the corresponding absorption spectra recorded under equivalent 

Q \ conditions. Ideally, one should also measure the induced terahertz (THz) response^i^iSiiSiiiii^. 

O i Based on our microscopic analysis, we determine criteria under which exciton population effects are observed and 

'>-J ' when electron- hole plasma contributions are dominant. We also discuss the Coulomb- induced nonlinearities observable 

• i-H in PL for elevated densities. For this purpose, one clearly has to perform a quantitative analysis where the microscopic 

^^ scattering, electron-hole plasma, and possible exciton populations have to be included microscopically. The developed 

'^ , analysis allows us to determine systematic criteria for a measured combination of absorption and PL to judge how 

■ - - ■ much of a given PL spectrum can be attributed to excitonic and/or electron-hole plasma emission. 

The discussion in this paper is organized as follows. We first present our general quantum theory for PL from 

coexisting electron-hole plasma end exciton populations. We show how this approach, with microscopic Coulomb 

scattering, can be formulated to compute the PL from quasi-stationary semiconductor populations consisting of a 

given exciton and electron-hole plasma population. We evaluate our theory numerically for a two-band QW model 

and present a detailed analysis of the exciton and electron-hole-plasma contibutions to the luminescence. We discuss 

different levels of simplifying approximations and emphasize the importance of the frequency and density dependent 

scattering. Finally, we summarize our results and use them to analyze experimental data. In the Appendices, we 

present technical details regarding the microscopic description of excitons and the density and frequency dependent 
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scattering calculations. 

II. MICROSCOPIC THEORY OF EXCITONIC LUMINESCENCE 

We formulate our microscopic theory for the incoherent emission of a planar quantum well. In our calculations, we 
concentrate on PL resulting from recombinations involving only two bands, i.e. the energetically lowest conduction 
and the heavy-hole band. Multi-band effects can in principle be included in the theory, but they complicate the 
notation and do not contribute to the essence of our analysis. In the numerical evaluations, we use typical GaAs 
QW material parameters and choose the QW confinement such that the 2s resonance is 6meV above the Is state as 
observed in recent experiments^ii^iiS^ For semiconductors, the microscopic description of carrier excitations can 
be obtained with the help of the fermionic operators ackn and a]^ j^ for conduction-band electrons and a^^kii and 

al j^ for valence-band electrons having momentum k|| in the QW plane. The quantized light field is described by 
bosonic operators Bq..q^ and B^ related to free space eigen modes with momentum q — {qn,q±) that is divided 
into the q|| component in the QW plane and q± perpendicular to the QW. The full system Hamiltonian of the 
interacting carrier-photon system is discussed e.g. in Refs. Ilmla and the relating microscopic PL dynamics is derived 
in Refs.l(Tl ll(lll8yi9ll2(ll2lll22L Here, we therfore present only a brief summary of the elements important for our analysis 
of excitonic PL. 

Since we want to compute PL from incoherent quasiparticle excitations, we concentrate on conditions where the 
coherent polarization and the coherent fields (Bq) vanish. We also assume that the system is homogeneously excited 
and quasi stationary. For such conditions, the single-particle distributions 

/k|| = («c,k||,cac,k||), /k|| = (ai-.kiia^ii), (1) 

for the electrons and holes can be taken as constant. 

The relevant two-particle correlations in the incoherent regime are excitonic and pure intraband carrier 
correlations pifi^ For the study of PL, we only need to evaluate the excitonic correlations defined as 

lll^klpkli _ A ; t t \ 

'^X — ^\'^c,k||'^t,,k]|"c,k]^+q||at),k||-q||/ 

= (ac,k||'3^t,,k]|"c,k]|+q||a«,k||-q||) - (ac,k||'3^t,,k]|"c,k]|+q||ai>,k||-q||)s, (2) 

where the subtracted single-particle contribution (a^aJ,ay'aA')s is obtained by using the usual Hartree-Fock approxi- 
mation {a\ala^'ay)s = {a\ax') {ala^i) - {a\ai,,) {alaxi) . 

The incoherent photon correlations and the microsocpic source for the PL are obtained from the photon-assisted 
electron-hole recombination and the photon-number-like correlations 

nk||,q||,gi = ^('^q||,9j_'^«,k|i"c,k||+q|,) 

== (^q||,9j_'^i;,k||"c,k||+q||) " (Sq|| ,j^ ) (fl^ j^.^^ ac,k|| -|-q|| ) , 

respectively. These are obtained by subtracting the classical factorizations from the full expectation values. Once 
these quantities are determined microscopically, we can evaluate the semiconductor PL and study the influence of 
electron- hole plasma and exciton populations. 
The exciton correlation dynamics follows from 

^^^^q„,k|,,k„ ^ ^q„,k|,,k,|^qii.k;i.k„^^q„,k;,k„ 
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where we have introduced the renormalized electron-hole pair energy 
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Furthermore, we have defined a collective photon assisted electron-hole recombination term 



(6) 

(7) 



where J\j = dvc\ -^u^^ determines the strength of the light-matter interaction via the dipole- matrix element d^c and 
the mode function Uq at the QW position. The square-root term is the so-called vacuum field amplitude containing 
the light frequency uj^. 

The single-particle source of Eq. Q originates from the Coulomb interaction and has the explicit form 



cqii.kjpkii ^ y 

J V'k,+q||-k|| 
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This expression describes the Coulomb-induced scattering contributions which develop the electron-hole plasma into 
a correlated state even when no correlations are present initially. 

Equation 10} also contains several Coulomb sums. Here, we have explicitly presented only those terms that contain 
the phase-space filling factor {I — f^ — f^) since they describe the attractive interaction between electrons and holes 
allowing them to become truly bound clcctron-holc pairs, i.e. incoherent excitons. These Coulomb sums are referred 
to as the main-sum contributions.^'^ The other Coulomb sums, denoted symbolically as -Drcst, consist of terms of the 
form {f^ — f'')J2^ c. Since {f^ — /") is much smaller than (1 — /'^ — /'*) for low enough densities, D^cst has a 
much smaller effect on the exciton correlations than the main-sum terms for conditions where bound excitons still 
exist. ^'^ Thus, we may apply the main-sum approximation where Drcst is omitted when we compute the excitonic PL. 
Besides the two-particle correlations, Eq. Q also contains three-particle correlation terms symbolized by T"^-^ ■^. The 
description of these terms is presented in Appendix B. 

The incoherent carrier quantities enter as sources to the semiconductor luminescence equation^ 
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is the renormalized kinetic energy. The first term in the second line of Eq. H1U|I shows that photon-assisted recom- 
bination processes can be seeded either by an electron-hole plasma population (/'^/'') or an exciton correlation term 
(cx)- Since these spontaneous emission sources do not depend on the photon frequency, both plasma and exciton 
correlation sources contribute to the PL in the same spectral range. 

The second term in Eq. (|10|l shows that the photon-assisted recombination processes for different momentum 
combinations are coupled via the Coulomb interaction. This renormalization determines the emission frequency 
similar to the classical field renormalization in the semicondcutor Bloch equationsii22i24i2^. Since this Coulomb 
correlation effect is not at all related to the specific spontaneous emission source, even the existence of bare electron 
and hole densities leads to emission of light at the bound exciton frequencies. When the QW is not inside a cavity, 
the stimulated feedback term, A{B'^ B^), in Eq. (|10|) is negligible. The explicit form of the three-particle scattering 
term, T^, is discussed in Appendix B. 



A. Excitons vs. correlated electron-hole plasma 

As a very broad definition, an cxciton is a quasi-particle that consist of a bound electron-hole pair. This interpre- 
tation is seemingly simple and straightforward. However, once one seriously considers the fact that an exciton has 
a fermionic substructure, as is done in Eq. |0J, one realizes that one can neither define a pure exciton operator nor 
a rigorous number operator for excitons in a fermionic many-body systemi2& As a working definition, we speak of 
an exciton in a many-particle system if at least one electron depends on a hole in a correlated manner, where the 
correlation is in the relative motion coordinate. The pair wavefunction is determined from the Wannier equation that 
is analogous to the generalized Schrodinger equation for the hydrogen. The bound state solutions define the bound 
excitons while the ionized states define ionized excitons. 

The Wannier equation in momentum space is equivalent to the Fourier transform of the Hydrogen equation, but can 
contain additional energy renormalizations and a phase-space filling factor if a finite amount of carriers is present tf 
As shown in Appendix A, it is even possible to incorporate carrier scattering in the Wannier equation and thus obtain 
excitonic quasi-particles including Coulomb scattering effects. In the presence of carriers, both the phase-space filling 
factor and the microscopic scattering make the corresponding eigenvalue problem non-Hermitian, such that a left- 
handed and a right-handed basis set (^i/.q,, ) must be introduced, see Appendix A. There, we show that the definition 
of true excitons can be mathematically formulated by using the transformations 

AA^.,.(q,|) = Y. <J>L„(k||)K„(k||))*4"''""''^''"^", (12) 

k||,kj| 

^ (<i>f,q,^(k|| -q,))*$^q,,(kj| +q,)AiV;,.,(q||), (13) 

where the excitonic center-of-mass momentum has been split symmetrically via 

qe = ■ q||, q/i = ■ q||. 14) 

nie + rrih nie + nih 

In particular, the diagonal elements, 

AiVA(q||) = AiVA,A(q||), (15) 

define the momentum distribution for excitons in state A. 

The physical relevance of these quantities can be verified by analyzing THz absorption which is shownifl*2i to probe 
exclusively ^ AA'^xlqii) under incoherent conditions. More precisely, each exciton state can be connected with a 
population 



q||,k||,k|| 



An^ = i^A7V,(q||) (16) 



qii 



and a ls-to-2p resonance in the linear THz absorption shows unambiguously that the Is state is initially populated. 
In full analogy to the optical transitions in atomic systems, the strength of the resonance is directly related to the 
population difference between the initial Anis and the final An2p states. Thus, THz spectroscopy shows that An^ 
defines the exciton density of a given state in the same way as one can count atoms occupying different atomic states 
by using optical spectroscopy. 

In general, exciton correlations can also exhibit nondiagonal X ^ v contributions in AA'^a.i^- As shown in Ref.^'i, 
these off-diagonal parts build up on a pico-sccond timescale and can be identified with a correlated plasma while the 
diagonal exciton populations are formed much more slowly. In the present publication we include microscopic Coulomb 
scattering beyond the constant dephasing and we include both diagonal and non-diagonal exciton correlations. 

B. Incoherent quasiparticle states 

As already argued irki^, it is essential to analyze the experimental PL together with absorption spectra obtained 
under identical conditions. Since many experiments measure PL in the normal direction, we consider this configuration 
here. Using our formalism, we establish a connection between the quasi-particle state and the measured spectra in 
order to determine the conditions and excitation parameters where excitonic populations significantly influence PL 
experiments. 



When a weak classical optical pulse, E{t), is applied to probe the incoherent excitation state of a semiconductor, it 
does not alter electron, hole, or exciton populations. The probe pulse interacts with the system, merely by inducing 
the optical (probe) polarization Pk™ = (a], k '^c,k«) in the system. Thus, we need to solve the polarization dynamics 
from the semiconductor Bloch equation o^ i 



d 



i-/k|| -/k|i Zl^s-'^ii^'' 



k;, 



i-fL-ft]dc,vE{t)~ir^.. 



(17) 



which contains two-particle Coulomb- and phonon-induced scattering terms in Fkn . Performing a Fourier transfor- 
mation of Eq. H17|) , we obtain 



(^w-ekii)-Pkii(w) 



l-/k 






- i-/k||-./k||]<.sM-*rk„(^), 

where we describe the Coulomb-induced scattering 

Tk,, {iv) = 7^„ {co)P^., ico) - E 7^°k;, (^)^k;, (■ 

k;, 



w 



(18) 



(19) 



via the diagonal and off-diagonal excitation-induced dephasing matrices 7^ and 7*^^ . Their detailed form is given in 
Appendix B. 

The linear absorption spectrum follows from 



a(Lu) — Im 



1 dk;, ^k;, (w) 

7Is 2^ 



eoS 



E{u) 



(20) 



If the semiconductor system is in thermodynamic equilibrium, the linear absorption, a{uj) and PL spectrum, PL{lu) 
are connected via™ the well-known thermal relation 






(c) 



a{uj) 



^{hLO-tJ.)/{kBT) _ I 



(21) 



where w is the angular frequency of the light, /i is the chemical potential of electron- hole pairs, and T is the lattice 
temperature. For dilute densities, the Bose-Einstein distribution can be replaced by a Boltzmann distribution. Clearly, 
Ip^{uj) provides a well-defined reference point under all excitation conditions. In fact, the relative contribution of 
the electron-hole plasma to the excitonic PL can be identified by quantifying how strongly the actually measured PL 
deviates from the thermodynamic limit. 

Under quasistationary conditions, the steady-state PL emitted in normal direction, qn = 0, is obtained from 
Eq. ©-|ini) as 



IpLiuJ = c|q|) = ^A(st_^^So,,J - Re 
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Here, we have omitted the stimulated contributions since they are insignificant for bare QW systems. As shown in 
Add. lEl the three-particle scattering terms can be cast into a form similar to the polarization scattering in Eq. H19|l . 
i.e. 



rpU 

-'k||,0,9i 



7^1 {^q^ )nk|| ,0,g^ - E 7^|°k; (w)nk' ,0,9^ ■ 



(24) 
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Thus, the semiconductor Bloch and luminescence equations exhibit the same generic structure. However, as the major 
difference, the semiconductor PL has as its incoherent source the combination of carrier densities f^f^ and exciton 
correlations while the optical polarization is driven by the coherent field E{uj). 

The incoherent quasi-particle correlations enter Eq. (|23l) via cx ■ To separate the exciton-population and correlated 
electron- hole plasma contributions, we introduce 



0,k||,k|| V^/^R „ A*^R /, /NA,r /„N . "'k;.k|| 



E('J'A:,„(k||)) <i>f:,„(kj,)AAr...(0)+c;;ii, (25) 

A 



where the steady-state form of the correlated plasma part is obtained from Eq. Q) 

9^ o,k;|,k|| _ o,k; k|| o,k;|,k|| 



+ Vk',-l 
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+ 1 - /k„ - /k^ E ^h -H siaLl - 1 - /k - /A, E ^'„ -■ - 



'Z^ [ 'k||,l||J Vasma * / . 'kj^ .1|| ^plasma- \'^^ I 

111 1|| 

Here, the main-sum approximation has been made and the spontaneous recombination term is omitted since it mostly 
influences the exciton population. We have also presented the three-particle scattering terms induced by the Coulomb 
scattering via 7/*" , , . 

^ 'k||,k|| 

The incoherent particle states are characterized by the respective single-particle carrier distributions /'^' '' and the 
exciton populations ATV^. Hence, under the assumed quasi-stationary conditions the PL can be solved for a given 
configuration by specifying Z*^' '' and AiV^. For low densities, and if we assume that the scattering can be approximated 
by a simple dephasing constant 7, Eq. l(^ reduces to 

A<r^(q||) = ^^^^%-E'J'A(k||)/k„+.Jk^-,.1>^(k||) (27) 

showing that this correlated plasma contribution does not contain densities since the diagonal terms v = \ vanish. 

C. Analytical approximations 

To obtain an analytical understanding of the fully numerical solutions of the complete equations, it is beneficial 
to introduce a generalized exciton basis in order to diagonalize the homogeneous parts of Eqs. (|18|l . %l'i\ and H26|l . 
Since these equations include also carrier-carrier scattering, energy renormalizations and phase-space filling factors, a 
suitable exciton basis clearly must go beyond the usual low-density Wannier wavefunctions. 

The basis set for solving Eqs. (|18|l and (|23|) is obtained from the solution of the generalized Wannier equation, 

EMk„k;,(^)¥'?(k'||,c.) = (ii;,(a,)-zr,.(c.))^^(k||,c.) (28) 

with the matrix M defined as 

Mk,„k;, (^) = [ek,i - 111, (^)]4;,,k„ - [l - /^„ - /k^] ^k;,-k|, + *7^°k;, {^). (29) 

Here, the scattering matrices 7^/*-"^(aj) are given in App.IbI Due to their occurence and because of the phase-space 
filling factor, Eq. (|28|l is clearly non-Hcrmitian such that we have to distinguish between left- and right-handed basis 
fmictions. The left-handed functions (pj(k||, cj) satisfy the equation 



E¥'^(k'||,^)Mk; k„(w) = (ii;.(a;)-zr.(u;))v5^(k||,w) (30) 



with the same eigenvalues as the right-handed ones. Since the scattering for the polarization and the photon-assisted 
polarizations is frequency dependent, also the resulting basis sets {ip^ (k|| , w)} and the complex eigenvalues Ei,{lo) — 
iVu{Lij) depend on lj. 

Introducing this generalized exciton basis in Eqs. H18I) and (I2U|I . we obtain the equivalent formulation 

a(w)(xdcv > -^r-i—\ — ^ -p / N A^A^), 31 

f,^ ^ ^ ^ ' 

with the driving term A^^fj_{uj) defined as 

= (^^Mi(i-r-/'>^H) (32) 

where we introduced the symbolic bra-ket notation for the inner product. 
In the limit of constant 7, the relation 

(l-/k„-/k„)^^(k||'^) = <(k||''^) (33) 

can be proven such that A^^^{uj) reduces to 

A,,^M =^(^^(k||,a;))*^;?(k||,c.)-5,,,. (34) 

and we recover the usual Elliott formula. However, in the present paper, we always solve the absorption with the full 
source term, Eq. (|S^ . 

Due to the structural similarity between the semiconductor Bloch Eq. H18|) and the luminescence Eq. (|23|l and the 
identical form of the microscopic scattering matrices, the photon-assisted polarizations nk,, g can be expanded using 
the same exciton basis functions. Inserting the solution of Eq. H23() into Eq. (|22|l . we obtain 

rPLf ^ |, ,2v^ ^^(r^0,a;)(y>g(r = 0,c^))* 

/ o) ex dcv 2^ ^r-i—\ % -r / \ ^'^A^)> (35) 

which is structurally identical to Eq. (|31|) and exhibits the same resonances and the same frequency dependent 
broadenings. The only difference lies in the PL source term 

5.,^(c.) = 5](^^(k||,c.))V/^„/£,^+5]4'''"''"]^^^(k||,u;). (36) 

k|| ^ kji ^ 

As we noted already in the discussion after Eq. p3|l . the PL source consists of two different contributions. The first 
one defines the Hartree-Fock exciton correlation 

S^ico) = {XlX,)ico)\^^ = ^ (^^(k||,c.))Vk„/k'„¥'^(k|h^)- (37) 

It contains both diagonal and non-diagonal contributions and is by itself not a consistent approximation for calculating 
PL spectra. Nevertheless, we will see later that it still gives a qualitative picture around the Is resonance. 

The second contribution looks exactly like the exciton transformation defined in Eq. (|13() . However, before proceed- 
ing further, we first look at the solution of the last remaining Eq. H2()|l . which determines the electron-hole correlations 
entering in the source of Eq. H36|l . This equation can also be simplified by introducing a suitable electron-hole pair 
basis, but the scattering for the cx-equation is different from the scattering entering in the previous two equations. 
In particular, the scattering matrices are not frequency dependent since we assume stationary excitonic correlations. 
Therefore, Eq. I|2()|) should be diagonalized with an exciton basis obtained by using the correct excitonic scattering 

matrices in the matrix in Eqs. (|29|l and (|28l) . We denote this second exciton basis set by <I>p (^ll)- Expressing the 
excitonic correlations cx in the PL source, Eq. (|36|l with this excitonic basis, we can rewrite the source term as 

AS.^c.) =J2{^^{u.)\^^,)A{xl,X,,){^^,\^j:iu;)) (38) 
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where the brackets denote either the expectation value or an inner product. Thus, the total luminescence source 
is given by the sum of Eqs. l|?7|) and (|SS|l . In the spirit of Eq. (P5)l . the correlated source AS^^fj_{uj) can describe a 
correlated plasma or true exciton populations. In general, we refer to plasma PL, if we include only the correlated 
plasma from Eq. (|26|l .to the PL source AS'i/^^(a;). True excitonic PL is obtained after we add additional excitons 
ANi, = A{XlX^) to the source term. 

In all our computations, we include the Hartree-Fock source as well as the plasma correlations including the correct 
microscopic scattering, unless stated otherwise. However, in order to gain some more physical insight to the meaning 
of Eq. H35|) and only for that purpose, not for the numerical evaluations, one may simplify the microscopic 
scattering by replacing it with a phenomenological decay constant 7. As a result, the exciton basis sets tp^ and $1, 
become identical and basis as well as eigenvalues become independent of co. Furthermore, we may apply the result 
from Eq. 1)2 7|l for the correlated electron-hole plasma. By inserting it into Eqs. I|35|) - (|36[l and using Eq. (|37|l . we 
obtain a simple Elliott-like formula for luminescence. 



E^iio) — hcj — i'-f 



with the source term given by 



kii 



^Elliott ^ y |^^(k||)|2/^ 4^ + A{xl,X^^o) = N^^iq = 0) + AN^q = 0). (40) 



We now see very clearly that excitonic PL can result either from excitons or directly from electron and hole densities. 
The form of Eq. H40|l is very useful when one envisions the differences of plasma and exciton PL. For example, Ref.llSl 
shows that plasma PL from Fermi-Dirac distributions produces an exponentially decaying tail for energies above the 
fundamental gap. This feature survives also when the scattering is fully included such that one can deduce rather 
accurately the temperature of electron and hole distributions directly from the PL spectra. 

III. NUMERICAL RESULTS 

The full microscopic form of the semiconductor absorption and luminescence equations shows that excitonic reso- 
nances can appear under quite different conditions. Clearly, Coulombic scattering can strongly broaden these reso- 
nances. We also notice that the scattering terms contain highly nonlinear contributions of the carrier densities. We 
therefore start by investigating these excitation-induced dephasing effects^SiSS on different excitonic resonances. 

A. Excitation-induced dephasing 

We first illustrate the influence of the microscopic scattering on the higher bound states by evaluating the fully 
microscopic absorption and PL spectra for various carrier densities. More specifically, we assume Fermi-Dirac distri- 
butions at 40 K for electrons and holes and take three very low values for the densities which are usually all considered 
to be in the "linear regime" . In addition to the microscopic Coulomb scattering, we have included a small background 
7 = 0.125 meV which results from the effects of acoustic phonon scattering and the purely radiative decay. Figure ^ 
presents absorption and PL spectra side by side showing that the Is-exciton resonance lies 8 meV below the band 
gap. For the densities used here, it is hardly changed and its width is basically determined by the background decay 
constant. But while the Is-exciton resonance is almost unchanged even for the largest density used here, the mag- 
nifications in the insets clearly demonstrate a very strong broadening of the 2s-exciton resonance. The absorption 
spectra show that the 2s resonance merges with the continuum already at around 10^ cm^^, a density which is usually 
considered very low in typical GaAs-type QW systems. The same trend is visible in the PL spectra where it is 
almost impossible to detect a clear 2s resonance, instead it merely marks the onset of the exponential decay into the 
continuum. This fact explains why the 2s resonance is often not cleary observable in experiments. 

The relation between excitonic resonances and nonlinearities can be investigated most easily if we define the separate 
contributions of Eqs. (|^ and l|S3|) via 

a{iu) = 5]a,(c.) (41) 

1/ 

I^^iLo) = ^/PL(^)(t^) (42) 



-8 



E 



4 - 



-■ — 1 — •- 


— T"^ 

0.4 - 


' 


- 


0.3 - 


■ 


1 
1 
1 




0.2 - 

-J 


,w~ 


■.A. 







■ 














0.3 


_ 




. 


6 


" 










d 


. 


0.2 


- 




fl 


JD 










1 1 


&_ 










V\ 


iL 3 


- 1 


0.1 




^\^ 


D- 


' 





. 





.J 


,i , , 


.i 1 



-8 



hv-Ee [meV] 



FIG. 1: Absorption (top) and PL spectra (bottom) for a carrier temperature of T = 40K for three different carrier densities 
of 71 = 5 X 10^ cm~^ (shaded area), n — 1 x 10^ cm~^ (sohd hne), and n 
magnification of the region around the 2s resonance. 



X 10^ cm ^ (dashed line) . The insets show the 
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(43) 
(44) 



Here, the individual contributions of the different exciton states v appear separately. However, since the eigenenergies 
and broadenings are frequency dependent, this simple form is in fact somewhat misleading. Due to the frequency 
dependence of the complex excitonic eigenvalues the generalized Lorentzian function can in reality deviate quite 
strongly from a simple Lorentzian — even more so since in general the oscillator strength in the numerator can be 
complex and thus mix real and imaginary parts of the energy denominator. 

In Fig. 13 we show the same three absorption spectra as in Fig. ^ as thick solid lines. But now, we also plot 
Oils, Oi2s and the rest separately. Since the excitonic states have a complicated co dependence due to the scattering 
matrices, we identified the Is and 2s resonances by choosing at each ui those two excitonic states which have the 
highest oscillator strength. Especially for higher densities, we find an entire continuum of excitonic states; some of 
them have a relatively small oscillator strength. But a clear distinction of the different states, as in the zero-density 
analytical Wannier formula, is no longer possible. 

Even for the lowest density, the 2s resonance is strongly non-Lorentzian and extends well into the continuum. 
Reversely, even the higher lying continuum states contribute significantly (almost a quarter) to the observed 2s peak. 
For the intermediate density, this effect becomes even more pronounced and for the highest density, it is no longer 
possible to determine a clear resonance. This becomes understandable if we look at the vertical lines in the figure 
which indicate the positions of the real parts of the eigenvalues. For low densities these real parts agree with the 
expected energies of the Is, 2s etc. states. However, we see that for elevated densities the microscopic scattering 
leads to resonances even between the Is and 2s energies. For the highest density, these resonances can extend down 
to the Is peak itself. 
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FIG. 2: Absorption spectra for a carrier temperature of T — 40 K for three different densities of n = 5 x 10* cm^"^ (top), 
n = 1 X lO" crn"^ (middle), and n = 6 x 10^ cm~^ (bottom). In addition to the full spectrum (solid line), the separate 
contribution for the Is resonance (dotted line), the 2s resonance (shaded area), and the rest (dashed line) is given. The vertical 
lines display the real parts of the excitonic eigenenergies of the generalized Wannier Eq. 12S 



Exactly the same analysis as before is performed for our PL spectra. For the intermediate density, the result is 
displayed in Fig. 13 All features discussed in Fig.|21are clearly visible also here. The excitonic positions are the same as 
in the absorption analysis, since the exciton basis for semiconductor Bloch and luminescence equations are identical. 
As we have already mentioned in the discussion of Fig.^a clear 2s resonance cannot be distinguished. It is composed 
of a whole cluster of resonances and has to be analyzed with a full microscopic theory. A zero-density Wannier basis 
is clearly insufficient to understand these intricate effects. 

Figure^presents an important part of the information necessary to understand the strong non-Lorcntzian lineshape 
features. Here, we show the frequency dependence of the broadening of the two bound states with the highest oscillator 
strength. For simplicity, we continue to call the second resonance a 2s resonance in analogy to the usual Wannier 
nomenclature. Our observation from the spectra are confirmed here; while the Is resonance is almost independent of 
frequency and dominated by the background decay constant 7bg = 0.125meV, the broadening of the second resonance 
is strongly frequency dependent. At the frequency value where the maxima of the corresponding resonances are 
observed (indicated by the vertical lines) the broadening of the second is already about five times as big as that of 
the first resonance. The fact that this broadening decays more quickly towards lower frequencies than towards higher 
frequencies explains the non-Lorentzian feature and the larger extension of the 2s resonance towards the continuum 
states. At the same time, the small Is broadening is consistent with the observed narrow Lorentzian line from 
Figs.HHSl 

This explanation covers only part of the non-Lorentzian lineshape. By looking at Eqs. H31() and H35|l we note that 
fundamentally the influence of scattering can be even more complicated. Since the presence of scattering generally 
leads to complex wave functions and correlations A^^^ and S^,^, even the oscillator strength can be complex and real 
and imaginary parts, i.e. broadened (5-function and principal value contributions in the complex Lorentzian functions 
can be mixed. 

Before turning our attention to more details of the PL spectra, we summarize the important finding of this section. 
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FIG. 3: PL spectrum and the different contributions for tlie same parameters as tfie middle frame of Fig.|21 
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FIG. 4: Imaginary part of excitonic eigenenergy of the generalized Wannier Eq. 1281 for the two bound states with the hightest 
(solid line) and the second-highest (dashed line) oscillator strength for a carrier temperature of T = 40 K and a carrier density 
of n = 1 X 10^ cm~^. The vertical lines indicate the respective real parts. 

We have noted that both in absorption and in PL the so-called "2s resonance" actually consists of a whole cluster of 
energetically close resonances. Except for ultra-low densities the usual nomenclature fails and the 2s resonance is not 
necessarily the second lowest bound state. Microscopic Coulomb scattering leads to the presence of a continuum of 
states, sometimes even relatively close to the Is-exciton resonance. 



B. Plasma luminescence spectra 



Since we have seen in the previous section that the analysis of the 2s resonance is a rather delicate issue, we start 
our detailed investigation of PL spectra by comparing different approximations and their influence around the Is 
resonance. This resonance is modified much less by Coulomb scattering than the 2s resonance such that we can 
anticipate it to be more robust under different approximations applied. 

In the present section, all PL spectra are computed from Eq. H35|l with the Hartree-Fock source, Eq. HH7(I . and 
the correlated plasma contribution obtained from Eq. l(^ . We do not add additional diagonal exciton populations 
until the next section where the ratio between excitonic and plasma luminescence is investigated. For comparison, 
we also present solutions of the simplified Elliott formula, Eq. (jSHJ, where the decay constant is chosen such as to 
match the width of the Is resonance of the full computation. Within the Elliott formula approach, there are three 
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different possibilities to deal with the four-point cx correlations. The best and most consistent way is the full solution 
of Eq. H2B|I with the scattering matrices replaced by a constant 7. The other two possibilities truncate already part 
of the four-point terms and can thus lead to inconsistencies at higher bound or continuum states. We include them 
here for comparison only and in order to relate previous results to the full theory. In the simplest case, we can neglect 
the Cx correlations completely. Then we obtain spectra in the same approximation as in our first publications^, 
referred to in the present paper as "Hartree-Fock" approximation. We can improve this approximation by solving 
the four-point correlations in a scattering approach as has been discussed in Sec. 4.4 of Ref. Ts. In that case, the 
four-point correlations are obtained as 



o,k;,k|, 



-S' 



0,k' kii 



,k' k„ 



(45) 



Z7 



with the definitions from Eqs. ([SJl and ||SJ|. As can be seen from the energy denominator, Coulomb resonances are 
not included in the description of the electron-hole correlations in that approximation. Hence, we have four different 
ways to calculate the Is PL: the fully consistent plasma result from Eq. H2t)|) with microscopic scattering matrices, 
the same with a constant 7 approximation (Elliott) , and then two approximations where either Eq. H26|l is truncated 
and solved without the Coulomb sums (scattering approximation) or where the four-point correlations are neglected 
completely (Hartree-Fock) . 
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FIG. 5: PL spectrum in the spectral vicinity of the Is resonance. The full result from Eq. I35II (shaded area) is compared to 
three different approximations of the Elliott formula with constant dephasing, Eq. 13911 . Here, the correlated PL source term 
is treated fully (solid line), in scattering approximation (dashed line), or set to zero (Hartree-Fock approximation, dark-shaded 
area). The constant-gamma results have been shifted by 40/ieV in order to collocate the peaks at the same frequency. 



Figure [3 presents the PL spectra around the Is resonance for these four different ways of computing the excitonic 
correlations. All approximations exhibit a clear peak at the Is resonance. A small shift of 40/ieV for the three 
constant-gamma calculations compensates small differences in the energy renormalizations with and without micro- 
scopic scattering. While the Hartree-Fock approximation gives only about one third of the full plasma PL, the other 
two approximations are within 30% of the full plasma value. Since we adjusted the decay constant, we can directly 
compare the peak heights in order to compare the quality of the different approximations. This is done in detail 
in FigEl for the two carrier temperatures of T = 20 K and T — 60 K. For both temperatures, the Eliott formula 
approaches the full plasma PL for very low densities, deviates most strongly at some intermediate densities, before 
getting better again in the limit of higher densities. The maximum deviation is temperature dependent and is less 
than 40% for 20 K and even smaller around 25% for 60 K. The Hartree-Fock approximation shows a similar functional 
form, but is in general too small by a factor of around three to four, with the precise value also depending on carrier 
density and temperature. The scattering approximation decays monotonously from low to high densities. For very low 
densities, it can even overestimate the correct value of the plasma luminescence as can be anticipated by extrapolating 
the curve to zero density. For high densities, it is too small by less than a factor of two. 
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FIG. 6: Relative peak height of Elliott formula (solid line), second-Born source (dashed line), and Hartree-Fock source term 
(dotted line) relative to the peak height of the full plasma result with microscopic scattering. Results for carrier temperature 
of T = 20 K (top) and T = 60 K (bottom) are shown. 
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FIG. 7: PL spectra for the same parameters as in Fig. |^ in semilogarithmic plot. The full plasma PL (solid line) and the 
Elliott-formula result are compared to the thermodynamic limit (shaded area). 



C. Excitonic plasma luminescence 



The results in Sec. IIII Bl show that the intensity of the Is-plasma PL can be predicted rather accurately using 
different levels of approximations. We next quantify the characteristics of the Is PL with respect to other resonances 
in the spectrum. Especially, we may use the thermodynamic relation (|21|l between a{u)) and PL(aj) to define how 
nonthermal the excitonic plasma PL is. The computed full (solid line), Elliot (dashed line), and thermodynamic 
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(shaded area) PL are presented in Fig. |7| Here, the carrier density has been taken as lO^cm"^ and the lattice 
temperature is 40 K. The shown thermodynamic PL is simply obtained by computing the microscopic absorption 
from Eq. (j2lj. Then, a((jj) is multiplied by a thermal factor according to Eq. 1)21(1 . The temperature is chosen such 
that the high-energy PL tails of the full and the thermodynamic computations agree. 

We observe now that event though the full and the Elliot PL agree at the Is resonance, the Elliot PL overestimates 
the high-frequency part of the spectrum since microscopic Coulomb scattering is not included. This is a clear indication 
that one needs to include the microscopic scattering to understand the quantitative features in the full spectral range 
of the excitonic PL. A comparison of the full and the thermodynamic PL indicates that the excitonic plasma PL is 
strongly nonthermal since it produces a much weaker Is resonance than thermodynamic arguments predict. We will 
see, later on, that one actually needs true Is-exciton populations to enhance the Is PL toward its thermodynamic 
value. The inherent nonequilibrium feature of the excitonic plasma luminescence has a simple explanation. It is 
intuitively clear that Is emission from plasma or correlated electron- hole plasma is more complicated than from 
exciton populations. This complication arises from the fact that an electron-hole plasma can emit a photon at 
excitonic frequencies, which are much below its average energy per particle, only by heating the remaining plasma. 
Due to this additional difficulty, the excitonic PL resulting from a pure plasma is generally weaker than that for the 
pure thermodynamic limit containing also exciton populations. 
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FIG. 8: "Phase-space" plot for /3-factor as function of carrier density and temperature. The temperature has been determined 
from the continuum tail of the emission. 

To quantify the level of nonthermal character in more detail, we define a (3 suppression factor. 



/3 



PL 



Is I full 



PLi 



s TD 



(46) 



first introduced by Schnabel et alm^. Essentially, the /3 factor tells us how strongly the actual Is PL, PLis|fuii, is 
suppressed with respect to its thermodynamic value PLi^Itd- Figure |S1 presents a contour plot of /3 as function of 
carrier density and temperature when only the fully mictroscopically included plasma luminescence is included. As a 
general tendency, the plasma PL becomes strongly nonthermal for low densities and temperatures. This is the regime 
which is favourable to sustain truly bound excitons. As temperature and/or density are increased, the excitonic 
plasma PL approaches the thermodynamic limit. 

As a general tendency, the (3 analysis has a rather weak sensitivity to the scattering model used or to low levels of 
disorder i2^ Thus, the /3 analysis is very much suitable to analyze the plasma nature of experimentally detected PL. In 
other words, the presented /3-phase diagram can be applied to diagnose the degree to which the measured PL in low- 
disorder GaAs-type samples is caused by plasma emission. On the basis of a quantitative analysis of simultaneously 
measured PL and absorption, Refsii^^ have shown that the experimental results for temperates above 30 K and 
densities above 5 x lO^cm"^ can be attributed to pure plasma emission. 
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D. Exciton-population effects 



Under suitable conditions a system of incoherent electron-hole quasi-particle excitations consists a mixture of 
correlated electron-hole plasma and exciton populations. Clearly, exciton populations may appear in all exciton 
states including the ionization continuum. In particular, the full thermodynamic equilibrium between plasma and 
excitons leads to a many-body problem^^i^ that produces highly nontrivial many-body states. However, several 
recent experimentai^*^i*^ demonstrate that excitonic PL is almost always far from the purely thermodynamic limit. 
Thus, it seems realistic to assume that the correlated electron-hole plasma coexists together with a fraction of excitons 
that is determined by the exciton formation kinetics and typically much below its thermodynamic value. Since the 
Is-exciton populations dominate at low temperatures, we include only them in our analysis of the exciton population 
contribution to the PL. 

As discussed in Sec. IIII CI the PL at the Is resonance can be strongly enhanced by the presence of true exciton 
populations. This can be seen directly from the luminescence Elliot formula H39|l . More precisely, a Is-exciton popu- 
lation basically adds luminescence at the Is frequency whose magnitude increases as a linear function of the exciton 
occupation AA^i^. In order to investigate this effect, we add Is exciton populations to the fully microscopically com- 
puted correlated plasma luminescence. We assume that the added Is excitons have a thermal momentum distribution 
such that AiVis(O), i.e. the 'bright part' of the exciton distribution that contributes to the luminescence in normal 
direction, is known. In practice, the population AAfis(O) is added to the source H25I) such that the fully microscopic 
luminescence is solved from Eqs. ()35|l for the coexisting exciton and plasma populations. 
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FIG. 9: PL spectra for a carrier density of n = 1 x 10^ cm~^ and a nominal carrier temperature of T = 20 K. The full plasma 
result (solid line) is compared to the thermodynamic curve (shaded area). Adding 1% excitons at a temperature of Tx = lOK 
(dashed line) drastically increases the Is PL. 



To illustrate the basic effect of exciton populations, we assume a low carrier density of lO^cm"^ and 20 K carrier 
temperature. Since the excitons can often reach a lower temperature than the carrier distributions, we choose a 10 K 
temperature for the excitons. Figure El shows PL without excitons (solid line), 1% excitons (dashed line), and the 
thermodynamic limit (shaded area). We notice that in all cases, the PL spectra are qualitatively similar, especially, 
the PL above the second resonance is virtually the same. Hence, the main feature of Is populations is that they 
increase the Is PL toward the thermodynamic limit. This confirms our previous assumption that exciton populations 
must be included to produce a thermodynamic PL. 

We investigate next how much the Is PL is changed by a 1% fraction of Is excitons as function of total carrier 
density. As in Fig. O we assume a 10 K exciton temperature. For each density, we define the ratio of the Is peak 
heights. 



Rx 



PLi, 



PLi 



s I full 



(47) 



between the case with excitons, PLi^jx, and without, PLi^lfuu, producing the microscopic plasma PL. The results 
are shown in Fig. ^| as function of carrier density for 20 K (solid line), 30K (dashed line), and 40 K (dotted line) 
carrier temperature. The exciton populations have strongest effect at low densities while Rx approaches unity for 
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FIG. 10: Ratio between Is PL with and without an added amount of 1% excitons. For an excitonic temperature of Tx ~ 10 K, 
the carrier temperature is T = 40K (solid line), T = 30K (dashed line), and T = 20K (dotted line). 



the elevated densities. This is in agreement with our earlier observation that high temperature plasma PL tends to 
approach the thermodynamic limit, i.e., excitons are then less important. However, we observe a strong increase of 
Rx for low densities and temperatures. According to Eq. H39|) . the low-density plasma PL scales with the product 
of the carrier occupations, f^f^, while the exciton population PL is directly proportional to the occupation AA^i^. 
When the carrier temperature is increased or the carrier density is decreased, the absolute magnitude of the plasma 
PL decreases nonlinearly. If the exciton occupation has a fixed temperature, as we have assumed in Fig. IIUI the 
relative strength of the population PL increases with respect to the plasma PL as the density is lowered or the carrier 
temperature is increased. Thus, the observed increase of Rx follows directly from the nonlinear dependence of the 
plasma PL for low densities. In particular, the exciton populations can contribute strongly under low temperature 
and low density conditions as verified in recent experiments>ii2^ 
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FIG. 11: Necessary exciton fraction to be added in order to obtain a ratio of 3/2 between the Is-resonance peaks with and 
without additional excitons. For an excitonic temperature of Tx = 10 K, the carrier temperature is T = 60K (solid line), 
T = 40K (dashed line), and T = 20K (dotted line). 



The effect of Is populations can also be investigated by evaluating how large an exciton fraction, 

^50% _ ^'^Islsoro 



neh 



(48) 



is needed to increase the excitonic Is luminescence by 50%. As in Fig. IIUI we assume that the Is-exciton distributions 
follow a thermal Bose-Einstein distribution at a temperature of 10 K. We assume carrier temperatures of 20 K (solid 
line), 40 K (dashed line), and 60 K (dotted line) and define x^° as function of density. The result is displayed in 
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Fig. ^] We observe again that low carrier densities are sensitive to the presence of exciton populations while their 
effect becomes less important for elevated temperatures and densities. The difference of low vs. high temperature can 
be understood on the basis of the same nonlinearity arguments as in Fig. E3| 

IV. DISCUSSION AND CONCLUSION 

Excitonic PL resulting from an electron-hole plasma has been observed in a variety of experiments performed with 
GaAs-type QW systems. For example, the ultrafast build up of the Is-luminescence resonance on a subpicosecond time 
scale REF is a clear indication of emission from a relaxing electron-hole plasma. The formation times of incoherent 
exciton poplulations are much longer. This has been carefully analyzed recently in experiments where absorption and 
PL spectra have been measured simultaneouslyii^*^ The microscopic analysis of these experiments not only verify 
that excitons are formed on a time scale of several hundred picoseconds but they also showed that appreciable exciton 
populations are obtained only for temperatures below 30 K and densities below 10^°cm~^. Otherwise, the excitonic 
PL could be explained by our plasma PL alone. 
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FIG. 12: Top: A comparison between the data points from Fig. 5 of Phys. Rev. B 71, 195313 (2005) and the thermodynamic 
result as obtained from analysing either peak heights (dashed) or areas (dot-dashed line) of Is and higher resonances, respec- 
tively. Our result without (dotted) and with (solid) added excitons obtained from the peak-height analysis has been added. 
Bottom: Necessary exciton fraction for an excitonic temperature of Tx = 10 K used for the top part of the figure. 



In order to check these interpretations, we analyze the experimental data published recently in Ref»i&. In Fig. 1121 
we have carefully read in all experimental data points from Fig. 5 in RefJ^ We compare these date to our full theory 
and several levels of approximation. We first compute the ls-to-2s ratio by integrating the thermodynamic result, 
Eq. (|21|l . obtained from our absorption spectrum for a density of 10^ cm~^. If we determine the ratio by integrating a 
spectral range of 8 meV around each resonance as has bene done in Ref. Il6l we obtain the dot-dashed line; if we plot 
the peak ratio instead, we get the dashed line. The area ratio lies well below almost all of the experimental data points. 
One can think of a large variety of possible explanations, which include radiative effects (the sample is measured in 
reflection geometry in front of a distributed-Bragg-mirror) or possible disorder effects which could potentially enhance 
the Is emission relative to the rest. In order to proceed with our analysis, we therefore use the peak ratio as a 
reference, which lies clearly above most of the experimental data points. The results of our luminescence theory 
for various temperatures is plotted in Fig. El before and after adding additional excitons to our plasma PL. The 
necessary fraction of excitons needed to obtain our fit is given in the bottom panel. Already the plasma theory nicely 
agrees with the experimental results for elevated temperatures above 40 K. And for lower temperatures we find that 
exciton populations do indeed contribute significantly. Even though we do not have access to the full experimental 
information, this comparison shows that the experimental results can be well explained with our theory, in contrast 
to the claims in Ref.[i^ For a detailed one-to-one analysis, the carrier density for each data point should be obtained 
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in the experiment as well. 

As an alternative experimental route, also THz experimentsii performed under similar conditions show unambigu- 
ously that Is-exciton populations are formed on a nanosecond time scale. Recently, also a combination of PL and 
THz absorption measurements has been reportedfi^ There it has been shown that a Is-exciton resonance can indeed 
be seen in the PL even when THz absorption does not exhibit a Is — 2p resonance indicating unambiguously that 
Is-exciton populations arc not present under these conditions. 

Both the strong non-thermal nature of PL spectra and the excitonic PL without a corresponding THz absorption 
signal can easily be explained if one includes the plasma PL in the analysis. While varying in importance, this 
contribution is fundamentally always present since it only requires that carrier distributions are excited. The strength 
of the present theory is to simultaneously treat effects from both the electron-hole plasma and the excitons by fully 
including the Fermionic substructure of the excitons. Only such a combined analysis without exclusive restriction 
to one of the possible emission sources can answer the question which source is important under given experimental 
conditona. 

To summarize, in the present paper we have discussed and evaluated our microscopic theory for calculating PL 
spectra of semiconductor QW systems and for analyzing exciton and plasma contributions to luminescence, respec- 
tively. We have particularly focused on the effects of microscopic Coulomb scattering and analyzed their influence on 
the different exciton resonances. While the PL Elliott formula, Eq. (|39|l . elucidates the underlying physics, a realistic 
theory-experiment comparison should always include microscopic scattering, in particular if the detailed structure 
around the band-edge is of special interest. A desirable extension of our theory would be to simultaneously include 
the exciton formation dynamics and the PL in a systematic way. Even though the theoretical framework exists, the 
numerical implementation of such a theory is very demanding and requires a lot of additional work and improved 
computer resources. 
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APPENDIX A: MATHEMATICAL DESCRIPTION OF EXCITONS 

Based on the four-point definition in Sec. lH Al we find a natural way to identify incoherent excitons once we evaluate 
the electron-hole pair-correlation function 

geh(rii) ^ (vl'l(r||)*[(0)*;,(0)*e(r||)), (Al) 

vl/e(r||)^-^^a,,k„e^''"-ii, vpt (r,,) = -^ ^ a^^^^^e'k,,-,,, (A2) 

I'll I'll 

where only the envelope part of the field operator ^A(r||) is included since the Wannier excitons depend on a length 
scale much longer than the lattice periodic part of the Bloch functions. In general, (7eh(r||) defines the conditional 
probability to find an electron at position r|| when the hole is positioned at the origin. Thus, properties of 5ch(r||) 
can identify what kind of correlations one can observe in the relative coordinate between the electrons and holes in 
the many-body system. 

The pair-correlation function can always be divided into single-particle as well as to correlated two-particle contri- 
butions and the homogeneous incoherent quasiparticle excitation conditions produce 



5oii(i"||) = neHh + A5oh(r||), (A3) 

k||,k]|,q|| 

We observe now that the single-particle factorization leads to a background contribution, rieTih, which states that 
probability of finding electron and hole simultaneously is proportional to the product of their respective densities. 
This contribution does not have a genuine r||-dependence. As a result, the true exciton effects are observed only in the 
correlated two-particle contributions Ag^hii^w) which contains exclusively the cx correlation. Thus, it is obvious now 
that only cx contributions can define and include the microscopic description of the true excitons in the many-body 
system. 
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The correlated Agoh(r||) can also be expressed in the exciton basis 






(A5) 



When A5oh(r||) exhibits a clear single-state r-dependence proportional to \tf^ (''n)!^' it is obvious that the ^_^ ^ sum 
is dominated by the element \ = v = \q. Since this case corresponds to having only Aq excitons, we may conclude 
that 



A"Ao = ^E^^^o,Ao(q||) (A6) 

defines the density of excitons in state Aq in the system. 



S 

qii 



APPENDIX B: MICROSCOPIC SCATTERING 

In order to get a realistic density dependence for the saturation of different excitonic resonances, it is important 
to include Coulomb scattering on a microscopic level. For completeness, we start by investigating the scattering for 
computation of absorption spectra, i.e. for the microscopic polarizations (Pfc) = (vlck). 

In the language of the cluster expansion, Coulomb scattering is provided by coupling of n-point expectation values 
to {n + 2)-point expectation values. In the case of the microscopic polarization, we find 

k' .q 

+ 51 ^9 \A{vlvl,Vy+qCk-q) - A{vl^^vlvk' +qCk)j . (Bl) 

k\q 

This equation is still exact, and the level of approximation is determined by the precision with which we calculate the 
four-point correlations. Their dynamics is determined by an equation of the general form 

z?i|a(4> =5[(2>]+T[(2),A(4>] + t/[A(6)], (B2) 

where the four-point correlations are driven by a source term 5'[(2)] containing only sums and products of two point 
terms (such as microscopic carrier occupations and polarizations), they are coupled to other four-point correlations in 
a complicated manner symbolized by T[(2), A(4)] and they are themselves coupled to yet higher order terms due to 
the hierarchy problem. The scattering approximation to this equation is obtained if all coupling between four-point 
terms and to six-point terms is neglected and only the kinetic energy contribution and the two-point source term is 
kept on the right-hand side. Formally, the scattering solutions can thus be obtained from 

in^A{A) ={e-i5)A{A)+S[{2)i (B3) 

Ot scatt 

where we have introduced a phenomenological dephasing constant 5 in addition to the sum of renormalized free-particle 
energies e. 

More specifically, we obtain 

d 

i^jrMAA'^^k'+qCk-q) = {eUq +£k+ Sk'+q " ^k' " i^) A{vlcl,Ck'+qCk-q) + Sk,k' .,q (B4) 



with the singlet source given by 

Sk.k',q = {Vq - Vq + k'-k) (/|.(1 - /|. + ,)(1 - Fk-q) + (1 " fk')fk'+q.fk-q) (Pk) 
-Vq (/|,(1 - /,V,)(1 - .4^) + (1 - rk').fk'+qfk) {Pk-,) 
+ Vq + k'-k (A^(l - .f^){l - fLq) + (1 - .fk')fk.fk-q) {Pk'+,)- (B5) 

Here, we have neglected source terms in third order in (Pj.) which is suitable for a low intensity probe pulse. In general, 
we want to calculate the susceptibility for a specific carrier density in a pump-probe type experiment relatively long 
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(up to nanoseconds) after the pump pulse. After the fast initial scattering, the carrier distributions are known to 
be changing very slowly and can thus be assumed constant. In that case, we do not have to apply the Markov and 
random-phase approximation as has often been done in the past. Instead, we see that the quantity Tk{Lo) which we 
introduced in the main part of the paper in Eq. (|17f) . is nothing but the Fourier transform of the right-hand side of 
Eq. (|B1I) . This quantity can be split into two parts, FJ: and F^, corresponding to the first and second line of Eq. (|B1|1 . 
respectively. 

Beginning with the computation of F^., we need the Fourier transform of Eq. IJB4|) which can be formally solved 
exactly via 



,t„t „ „ \c, ,^ _ Sk,k',q{^) 

'-k-q ^k ^k'+q 

In a similar manner, we obtain 



A{v,c„Ck,^,Ck-q)iu) = ^_ , ^ ^ ^^^ . (B6) 



A/ t t \/ \ Sk+q,k',q{'^) ,-n„. 

' ?iW - £fe - 4+q - ^k'+q + ^k' + ^S 

With these ingredients, we can calculate FJ: by summing over Eqs. ljB6p and (|B7p . using the definition of Eq. (|B5|I . 

Vk-q 



n= T. 



k'.,q ^^ ~ ^1 ~^k~ £fc' + ^I'+q-k + ^^ 



[{Vk^, - Vk>-k) [K'+q-kil - fk')a ~ fq) + (1 - fk'+q-k)fk'fq) (Pk) 
-Vk-q {fk'+q-ki^ fk')i^ Ik) + (1 - fk'+q-k)fk'fk) {Pq) 

+ Vk'-k (/feV,-fe(l - ft)a - fq) + (1 - fk'+q-k)fkfq) {Pk')} 

_v ^^ - 



^g ^^ - el -^q- ^l' + el'+fc-g + ^5 



[{V,.k - %'-,) (AVfc-,(l - A-')(l - m + (1 - fk'+k-q)fk'fk) (p,) 

-Vk-q {rk'+k-qil - /l')(l - /,') + (1 - fk' + k-q)fk'fq) (Pk) 

+ Vk'-q {rk' + k-q^ - fq)a " fk) + (1 " fk' +k-q) fq fk) {Pk')} (B8) 

The other quantity F^ can be obtained by symmmetry considerations. By exchaning each electron and hole index, 
we find 



k',q ^^ " ^« - ^1 - ^k' + ^fe'+«-fe + '^ 



{(Vk-q - Vk'^k) {fi>+q-k{l - /'')(1 - fq) + (1 " fk'+q-k)fk'fq) (Pk) 
-Vk-q {f^>+q-ki^ fk')i^ fk) + (1 - fk'+q-k)fk'fk) (Pq) 

+ Vk'-k ift'+q-ki^ - /fe)(l - fq) + (1 - fk'+q-k)fkfq) {Pk')} 



{ 



k',q ^ ~ ^k - el - et + et.+fe_, + ^6 



(Vfc - Vk'-q) {f^^+k-qa - /fc')(l - fk) + (1 - fk'+k-q)fk'fi) (P,) 
-Vk-q (/fe^ + fe-,(l - ft')i^ f'q) + (1 - f]:.+k-q)fk'fq) {Pk) 

+ Vk'-q {fb + k-q{l fq){l fk) + (1 " fk' + k-q)fqfk) (Pk')} ■ (B9) 

The previous two equations can be rewritten in terms of microscopic scattering matrices 7jJ'(w) and 7^?*}?/ (w) which 
couple microscopic polarizations with different fc-indices with one another. The scattering contribution to the polar- 
ization dynamics is then given by 

F, = -^^fiu;){Pk)+^Y,^^^,{cJ){Pk-), (BIO) 
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where the scattering matrices can be defined by comparison to Eq. IJB8|) and ljB9|) . While the coupling of the different 
components of the microscopic polarization leads to a redistribution among the (Pj.) and thus to dephasing, we can 
also see from Eqs. (|B8 | - (|B9)| that 



dt 

k 



scatt 



In terms of the scattering matrices, this relation is guaranteed due to their property 

7?M=5]7.^!iH. (B12) 

k' 

We want to emphasize that for stationary densities and a weak probe pulse it is natural to solve the equations using 
the Fourier transformation. In that way, we obtain frequency dependent scattering matrices where lo enters in the 
energy denominator. This is different from the usual procedure of using Markov and random-phase approximation in 
the time-domain. We can recover the old Markov result if for each term we replace the true frequency w by the single 
particle frequency us which describes the main evolution of the respective source. 

1. Scattering for photon-assisted polarizations 

For photon-assisted polarizations, the proceeding is quite similar to the previous subsection. In this case, the 
photon-assisted two-point quantities couple to photon-assisted four point quantities in the form 



k' .q 



+ Y.^i {^i^H'^W'+qCk-q) - A{b\vl^/^,Vk:+,Ck)) . (B13) 

k' .q 

As before, one has to set up equations of motion for these photon-assisted four-point terms which couple to yet higher 
correlations. In order to obtain a scattering approximation of the photon- assisted 4-point terms, we only keep the 
renormalized kinetic energies and the main sources. These sources are similar as in Eq. (jB5|l . but due to the presence 
of the photon operator many more factorization possibilities can occur. Under fully incoherent conditions, however, 
where we require expectation values of coherent quantities such as (6^), (Pfc) or {h^c^c) to vanish, the only remaining 
terms have exactly the form of Eq. IjBSp where each polarization has to be replaced by the corresponding photon 
assisted polarization. 

In the case of PL, we expect a steady-state solution since the source term of PL changes adiabatically and can 
be assumed constant. Nevertheless, the photon operator in each term provides a frequency dependence via the free 
rotation of the photon operator via 

ih^A{b\d) = -hu:,A{h\d) (B14) 

ut cm 

for each carrier operator O. When we set the left-hand side to zero, we find the same frequency dependence as for 
the Fourier transformation of the polarization equation. 

Consequently, the scattering matrices for photon-assisted polarizations are identical to those for the microscopic 
polarization and the scattering denoted by T^^q in Eq. H1U|I in the main text can be expressed as 

Tfc,, = -z7fe°K)(6tA.) +*E7°fc'K)(&Jft'>, (B15) 



k' 



in full analogy to Eq. HB10|) . 

2. Scattering for excitonic correlations 

For the calculation of PL we need to know the excitonic correlations within the radiative cone, i.e. with infinitesimal 
center of mass momentum. In the main text they are called 

A{civl,ck'vu) = A{clv^vl,ck') = A{PlPk'). (B16) 
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At this point, we can briefly address different approximations for the excitonic correlations. If we completely neglect 
excitonic correlations in the PL source term of Eq. (|35|l . we recover the old Hartree-Fock resullj^. If instead we set 
up the equations of motion for the excitonic correlations in the form of Eq. ljB2p and treat them on a scattering level 
by neglecting all coupling to four- and six-point terms, we can formally solve for cx and express the correlations 
as functional of carrier densities. In the main text this is refered to as scattering approximation of the excitonic 
source terms. Since photon number operators and photon-assisted correlations themselves are four-point operators, 
both approximations are fundamentally inconsistent since they neglect other four-point terms. This can be seen 
as negative PL at the continuum (scattering approximation) or even at higher bound states (Hartree-Fock), but 
nevertheless a qualitatively correct picture is obtained around the Is-resonance. 

In a consistent approximation, all four-point terms have to be treated on the same footing. Thus, the two-point 
sources and the four-point terms of Eq. HB2|I have to be solved and approximations are only introduced on the next 
level, for the six-point quantities describing scattering of excitonic correlations. Thus, we have to set up equations of 
motion for six-point quantities of the form 

in-^A{6) ={e-i5)A{6)+S[{2),A{4)l (B17) 

Ot scatt 

where we keep the free-particle energy and the source term consisting of products of two and four-point expectation 
values. Its precise form results from factorizing eight-point terms as e.g. {c\vyCy,c\ci+qiCk"+q~-q'Ck'Vk-q) ■ It turns out 
that all factorized sources contain at least one four-point correlation. Due to the indistinguishability of the quantized 
carriers, the general factorization procedure analogous to the previous two subsections leads to a huge variety of 
different triplet scattering terms. In particular, carrier operators belonging to one polarization operator in Eq. (|B16|) 
can be split into two different expecation values. The full Coulomb scattering including all those source terms includes 
all quantum effects due to the Fermionic character of the excitons as well as possible momentum scattering between 
different excitons. This full scattering scenario is currently not feasible in two dimensions. Therefore, we make 
the assumption that scattering acts only within one of the polarization operators of Eq. (jB16|) . In analogy to the 
previous section of photon-assisted polarizations, we assume that the second polarization operator remains unchanged 
and is not split into different expectation values. Under this assumption, the derivation of the scattering proceeds 
analogously as before, the number of scattering terms is also identical, and the only major difference is the lack of a 
frequency dependence of excitonic correlations. Thus, we obtain exciton scattering of the form 

"■^dl^x = -^[bk\ +lk')Cx 

I -X^/r-ODl* 0,k'.k" , ~D 0.fe",fc\ /iDTO\ 

+*Z^(bfc.fc"J Cx +Jk',k"Cx )■ (B18) 

fc" 
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